Avalanches, Breathers and Flow Reversal in a Continuous Lorenz-96 Model 
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For the discrete model suggested by Lorenz in 1996 a one-dimensional long wave approxima- 
tion with nonlinear excitation and diffusion is derived. The model is energy conserving but non- 
Hamiltonian. In a low order truncation weak external forcing of the zonal mean flow induces 
avalanche-like breather solutions which cause reversal of the mean flow by a wave-mean flow inter- 
action. The mechanism is an outburst-recharge process similar to avalanches in a sand pile model. 
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In 1996 Lorenz suggested a nonlinear chaotic model for 
an unspecified observable with next and second nearest 
neighbor couplings on grid points along a latitude circle 
[l[ . Due to its scalability the model is a versatile tool in 
statistical mechanics [Jlsl and meteorology The 
nonlinear terms have a quadratic conservation law and 
satisfy Liouville's Theorem. For strong forcing the model 
shows intermittency [9*]. 

The Lorenz-96 equations for the variable Xi are a sur- 
rogate for nonlinear advection in a periodic domain 



—Xi — Xi^i[XiJ^i — Xi^2\ — iXi -\- Fi 
at 



(1) 



7 characterizes linear friction (7 = 1 in 1]) and is a 
forcing. 

In this Letter a continuous long wave approximation of 
the Lorenz-96 model is derived. A surprising finding is 
that the nonlinear terms in the Taylor expansion are as- 
sociated with generic dynamic operators. Furthermore, 
the dynamics in a truncated version reveals avalanches, 
breather-like excitations and flow reversals, which mimic 
various physical processes in complex systems in a sim- 
plistic way. 

Lorenz [10] has analysed the linear stability of the 
mean m of Xi in ([T]) and found that long waves with 
wave numbers k < 2tt/3 are unstable for a positive mean 

TO. 

For = F = the equations ([T]) are conservative with 
the conservation law, Hx = ^f^^i^f, denoted as en- 
ergy in the following. The dynamics in the state space of 
the Xi is non-divergent thus satisfying Liouville's Theo- 
rem, dxjdx, = 0. 

The dynamics of an observable function Q{X) is given 

by 

Qt^{Q,Hx} (2) 

with an anti-symmetric bracket 

{A, B} = d,AJ,jdjB = -{B, A} (3) 

and the antisymmetric matrix 

Jij = Xi^iSj^i+i — Xj^i6ij+i (4) 

Energy Hx is conserved due to the anti-symmetry of the 
bracket. 



The conservative terms of the Lorenz-96 equations 
([1]) are obtained for Q = Xi. The equations are non- 
Hamiltonian llj since the Jacobi identity 



dJjk 



dJk 



ki 



dX, 







(5) 



is not satisfied. 

A continuous approximation is derived for a smooth 
dependency of Xi on the spatial coordinate x = ih'uv the 
limit /i — >■ 0. The variable Xi is replaced by a continuous 
function u{x, t) which is interpreted as velocity in the 
following. We use the infinitesimal shift operators 



E 

fe=0 

to write the bracket ^ as 



{±hd, 
kl 
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with 



(6) 



(7) 



(8) 



where {L^u) is a multiplication operator. The bracket is 
anti-symmetric since the adjoint is L*^ = L_. 

By taking n-th order truncations of the operators L±, 
we can find a hierarchy of truncated anti-symmetric op- 
erators 



where 



L±,n = ■ 



k=0 



kl 



The total energy for the velocity u{x, t) is 

1 



H = 



'dx 



(9) 



(10) 



(11) 



To each of these truncated operators corresponds a con- 
tinuous Lorenz-96 model 



Ut = {u,'H}r 



(12) 
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where the indices indicate the operator Jnm (as in ([T]) 
periodic boundary conditions are assumed). 

The expansion of the nonhnear terms in (JlJ up to order 
0{h'^) yields for the rescaled coordinate x' — ~x/3 (the 
prime is dropped below) 



Ut 



+ -UUr 



+ f 



(13) 



with an advection and further nonlinear terms which are 
due to the noncentered definition of the interaction in 

©■ 

The nonlinear terms are associated with antisymmetric 
evolution operators 



0{h) : Ji - -^(ud^ + d^u) 
Thus the evolution equation (1131) can be written as 



Ut = J 



5u ' 



J = Jl+j2 



(14) 
(15) 

(16) 



Note that the 0{h^) expansion in (fT2|) is represented by 
a third operator = —{l/18){uxxdx + dx^xx)', here we 
restrict to the 0{h?) expansion ([T3| . 

The evolution equation (I13|) has a conservation law 
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(17) 



(18) 



with the conserved current (j) which leads to the conser- 
vation of total energy ([TH) . Further conservation laws 
could not be found. In particular momentum given as 
the mean flow 



Thus, mean flows with [/ > (f/ < 0) are unstable (sta- 
ble) as in the discrete system ([T]) analysed in [lo| . 

The equations ([20l [22l) represent a coupling between 
perturbations and the mean flow. A forcing drives the 
mean flow towards positive values which allow the growth 
of perturbations. When the perturbations are sufficiently 
intense they reduce the flow to negative values causing a 
decay of their intensities. 

The nonlinear energy cycle represented by the ex- 
change between zonal flow and wave energy in ()20|) and 
is analysed in a spectral model for the unstable long 
waves by Fourier expansion in a periodic domain [— tt, tt] 



N 



a„ sm[nx) -f &„ cos(nx) 



ri=0 



Here we restrict to the low order system N = 2. 

ill = bobi + i&oai + 7; («i«2 + 61^2) 
D 2 

+ ^ (ai&2 - bia2) 

bi = -boai + ^bobi - \ (0102 + 6162) 
6 4 

+ ^ (ai&2 - bia2) 

2 1 

0.2 ^ 260&2 + -^boa2 + 2 (^1 + ^1^1 ~ «i) 

1 2 

62 = -2&oa2 - aibi + - (5? - aj) + -6062 



(23) 

(24) 

(25) 

(26) 
(27) 
(28) 



The mean flow is U = bg which is subject to a con- 
stant forcing / in the numerical experiments (|24p . The 
truncated system conserves energy 



Htnt — Hn 



H1+H2 



(29) 



U — lu) — / udx 



(19) 



is not constant. 

In the following we consider a constant and positive 
forcing / (note that the system is not dissipative). In the 
presence of perturbations v to the mean flow, u = U + v, 
the mean flow energy H — U"^ /2 changes according to 



The perturbation energy 



grows for positive U 



(20) 



(21) 



(22) 



hi Hi^\{a\ + b\) 



H2 = -(al + bl) 



The Liouville Theorem is not satisfied 




(30) 



(31) 



The expansion and contraction of the state space volume 
is controlled by the sign of the mean flow. 

Numerical experiments reveal a flow reversal mecha- 
nism and vanishing long term means of mean flow and 
wave number amplitudes, hence the Liouville Theorem 
(|3T1) is satisfied in the mean. 

(i) Weak forcing with / = 0.1 in the = 1 trunca- 
tion reveals periodic flow reversals (Fig. [Ij. The system 
starts with randomly chosen amplitudes. A mean flow 
increases gradually to positive values where it becomes 
unstable due to the excited waves, denoted as breathers 
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FIG. 1. Weak forcing / = 0.1: (a) Amplitudes 6o,ii,fei for 
= 1, intervals ~ 110, (b) amplitudes during a flow reversal 
with Eq. H33|) for feo (dashed), (c) energies Ho, Hi and Eq. 
(|34|l for Hi (dashed), (d) amplitudes 6o,a2,b2 for N = 2, 
(ai,6i vanish), and (e) energies {Hi vanishes). 



and the perturbation energy is 



cosh (at) 



(34) 



where a is related to the total energy H = 18a^. Hi 
attains its maximum during flow reversals when U = 0. 
These approximations are compared to the forced simu- 
lation in Fig. [T1d,c centered at a single flow reversal. 

In the presence of forcing / and for a small wave energy 
Hi the mean flow bo grows linearly in time, bo{t) ss ft, up 
to a value bo^max- This defines an interarrival time scale 
of flow reversals, r — 26o, mas//- In this range the wave 
energy evolves rapidly according to Hi{t) ^ exp(/t^/6). 
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FIG. 2. Intermediate forcing / = 1: Energy distributions 
(a) N =1, intervals ^ 13, (b) iV = 2. 



for 



The described flow reversal mechanism is retained for 
viscous dissipation represented by a linear damping of 
the wave amplitudes ai and bi. 

For the N = 2 truncation with all modes 6o, ai, 6i, 02 
and &2 flow reversals occur on a time scale roughly twice 
as for = 1 (Fig. [TJi). Due to the weak forcing the 
energy cascades to mode 2 with negligible amplitudes 
ai,bi and energy Hi (Fig. [Hi, e). Neglecting the modes 
1, the energy cycle for interactions among 60,12 and 62 
is 



in the following. These breathers drive a rapid flow rever- 
sal towards a negative flow which initiates their collapse. 
The process is energy conserving on short time scales. 
The total energy increases (decreases) when the mean 
flow is positive (negative). 

For A^ = 1 with the amplitudes 60,01,61 the energy 
cycle is for / = (compare ^ ^) 



dtHo = -^boHi, dtHi = hoHi (32) 



which is controlled by the mean flow. The solution for 
the mean flow is for 60 (0) = 



dtHo 



-boH2, dtH2 



:boH2 



(35) 



60 = — 6atanh(ai) 



(33) 



This corresponds to a rescaling of the Hq — Hi cycle ((32)) 
hy t — 2t for time and 60 — 26o etc. for the amplitudes, 
hence the energies quadruple. 

(ii) For intermediate forcing with / = 1 the time scale 
between flow reversals decreases by an order of magni- 
tude in the A' = 1 truncation (see Fig. [5^). Thus the 
intervals r approach the duration of individual breathers. 
For the complete set of modes in A^ = 2 (Fig. ^jp) the 
system is weakly nonlinear with a mixing of frequencies, 
uj/2,uj,3uj/2 and 2uj, where uj ~ 27r/T is defined by the 
interarrival times of the flow reversals ^12j. The lowest 
frequency determines the amplitude modulation. 

(iii) For strong forcing, / = 10, the flow reversals in 
the A^ = 1 truncation are regular (Fig. [3^) with intervals 
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FIG. 3. Strong forcing / = 
N = 1, intervals ^ 3, (bj iV : 



10: energy distributions for (a) 
= 2. 



decreased by an order of magnitude relative to / = 1. 
The dominant part of energy is accumulated in waves. In 
the N = 2 truncation the dynamics becomes intermittent 
as in the regime behaviour detected by Lorenz Q in the 
discrete equations (P). The events loose their identities 
and the systen becomes strongly nonlinear. 

In Summary, a continuous dynamical equation derived 
from the Lorenz-96 model is able to mimic several types 
of complex processes observed in geophysics, geophysical 
fluid dynamics, and solid state physics: 

(i) Avalanche processes excited by continuous driving 
as in the sand pile model of Bak et al. [l^; see also 
the recent observation of quasi-periodic events in crystal 
plasticity subject to external stress [T3|. A common char- 
acteristic property is the weakness of the external forcing 
which is necessary to cause avalanches. In the present 
model the flow is driven by a constant forcig towards a 



state where mean flow and wave energy interact. The in- 
tervals between the flow reversals are approximately pro- 
portional to the inverse of the forcing intensity, ~ 1//. 

(ii) The Quasi-Biennial Oscillation (QBO, [lal), a flow 
reversal in the tropical stratosphere driven by two dif- 
ferent types of upward propagating gravity waves. A 
common aspect is that the driving of the mean flow by 
waves occurs only for a particular sign of the mean flow. 
Although the QBO is considered to be explained dynam- 
ically the simulation in present-day weather and climate 
models necessitates careful sub-scale parameterizations 
or high resolution models [l^. The present model is 
clearly an oversimplification but can be considered as a 
toy model for this phenomenon. 

(iii) Rogue waves (also termed freak or monster waves) 
at the ocean surface are simulated mainly by the nonlin- 
ear Schroedinger equation (e.g. [13, EO; a Lagrangian 
analysis has been published recently [l^. The breather 
solutions found in the present model show characteris- 
tics like the rapid evolution and the high intensity in an 
almost quiescent medium. 

Due to the flow reversals the total energy of the non- 
dissipative system remains finite for a constant forcing. 
The long term mean of the mean flow vanishes and the 
Liouville Theorem ([5Tt is satisfied in the mean. The flow 
reversals are insensitive to viscous dissipation. 
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